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Abstract 

Light diffraction through a subwavelength aperture located at the apex of a metahic screen with 
conical geometry is investigated theoretically. A method based on a multipole field expansion is 
developed to solve Maxwell's equations analytically using boundary conditions adapted both for the 
conical geometry and for the finite conductivity of a real metal. The topological properties of the 
diffracted field are discussed in detail and compared to those of the field diffracted through a small 
aperture in a flat screen, i. e. the Bethe problem. The model is applied to coated, conically tapered 
optical fiber tips that are used in Near-Field Scanning Optical Microscopy. It is demonstrated that 
such tips behave over a large portion of space like a simple combination of two effective dipoles 
located in the apex plane (an electric dipole and a magnetic dipole parallel to the incident fields 
at the apex) whose exact expressions are determined. However, the large "backward" emission in 
the P plane - a salient experimental fact that remained unexplained so far - is recovered in our 
analysis which goes beyond the two-dipole approximation. 
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I. DIFFRACTION AND NEAR-FIELD OPTICAL MICROSCOPY 



Diffraction and scattering are among the most difficult phenomena encountered in 
optics. Eigenvalue problems with complex boundary conditions have to be solved in order 
to account for the finite wavelength of light in the interactions between optical waves and 
matter. Exact solutions are extremely rare and limited to simple, idealized diffracting 
objects. 

In general, approximation methods are used for diffraction problems involving apertures 
or obstacles of arbitrary geometry. For example, the interaction of light with a screen 
in the limit where the wavelength is short compared to the dimensions of the obstacles 
was already treated a long time ago by Huygens, Fresnel, and Kirchhoff (among others). 
This case corresponds to small deviations from geometrical optics. However, it is totally 
irrelevant for phenomena implying strong near-field patterns which occur for wavelengths 
comparable to the dimensions of the diffracting object. In these cases, the interaction of 
light with matter cannot be considered as a small perturbation, and the related boundary 
problem must be solved exactly. 



This is the situation encountered in Near-Field Scanning Optical Microscopy (NSOM) 
in which strong interactions exist between a nanosource and nearby nano-objects iH^. In 
NSOM, the specimen is usually illuminated through a subwavelength aperture located at 
the apex of a coated, tapered fiber tip. Such an optical tip can be mimicked by a truncated 
metallic cone, i. e. a very complex geometry. The aim of the present paper is to study 
this diffraction geometry in order to elucidate the properties of actual tapered optical tips. 
Although a rigorous solution of the problem will not be given, we will make use of an approx- 
imation method based on a field expansion in "quasi" spherical multipoles which is valid in 
the far-field. Maxwell's equations are solved analytically within this scheme, which allows 
us to describe in detail the topological properties of the diffracted field. In particular, we 
will show that the dipole behavior anticipated earlier js, 6| can be formally justified and that 
it remains valid over a large portion of space. The predictions of our analysis are found to 
be in a fair agreement with available far-field emission patterns of actual conical optical tips. 



The paper is organized as follows. In section II, the problem of diffraction through an 
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aperture in a metallic screen is introduced in general terms. In section III, a method to solve 
electromagnetic problems with conical geometry is developed. The radiation pattern of a 
conical optical probe is described in section IV for the case of a tip coated with a perfect 
metal, which is extended to real metals in section V. Section VI compares the theoretical 
results with experimental emission patters of metal coated, optical fiber tips. A summary 
is given in section VII. 



II. INTRODUCTION TO THE PROBLEM OF DIFFRACTION BY SMALL 
APERTURES : THE BETHE-BOUWKAMP SOLUTION 



In classical optics, the diffraction of light by an aperture in a screen is studied within the 
Kirchhoff approximation [3, ^. The use of this scalar method is justified in numerous cases 
where the polarization of light can be ignored and where the aperture size is large compared 
to the wavelength of the incident radiation. In order to apply the Huygens-Fresnel principle 
to light diffraction by small apertures, Smythejs, [lo| has developed a vectorial Kirchhoff 
integral formula based on the Green function method. The formalism is very general but 
- like Kirchhoff's scalar theory - its practical application is restricted to long wavelengths 
compared to the geometrical parameters of the aperture. The first-order solution consists in 
neglecting the influence of boundary effects on the light in the aperture zone. This method. 
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12| approximation 



which may be viewed as the electromagnetic counterpart of the Born 
in quantum mechanics, is adapted to the quasi geometrical regime of Maxwell's equations. 
However, it is not adapted to long wavelength radiation. 



Bethe 



13| (and after him Bouwkamp [l4|) have developed a rigorous solution for the case of 



a small circular aperture in a perfect metallic plane by solving an electrostatic and a magne- 
tostatic eigenvalue problem for the static near-field existing in the vicinity of the aperture. 
This solution, which generalizes previous work by Rayleigh 15|], can be easily recovered by 



using oblate spheroidal coordinates and harmonical wavefunctions. The Helmoltz equation 
[V^ + k'^]tp = 0, which reduces to the Laplace equation for small hole sizes V'^ip = 0, has 
the solution ip {(,,1], (p) = ■ PJ^ {irj) .e^^"^'^, where denote Legendre wavefunctions 

indexed by two integers l,m. C,, ">],(/) are the oblate spheroidal coordinates. The Bethe so- 
lution exhibits a complicated behavior in the vicinity of the aperture. Nevertheless, the 
Bethe solution reduces to a multipole expansion in the near- field domain limited by r /A ~ 1 
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(which corresponds to r — )■ +00 for the static Laplace equation where A — )■ +00). In the 
particular case of an incident plane wave, it reduces to a electromagnetic dipole field (see 
Fig. d]). The corresponding effective dipoles P^g and M^g depend on the aperture radius 
and on the incident electromagnetic field Eq, Bq as given by the formula P^g = =F^a'^Eo,_L, 
Mgg = ±Jj:a^Bo,|j. The signs in front of the expressions refer to the "incident + reflected" 
domain and to the "transmitted" domain, respectively. It is interesting to note that the 
factor of 2 between M^g and Pgff appears in all similar situations as we will see throughout 
this paper. The physical signification of the Bethe dipoles can be understood using the 



Clausius Mossoti 



16| formula. It describes the polarization produced by a locally constant 



electromagnetic field Eg, Bq in a small dielectric sphere of constant permittivities e, /i which 
is immersed in a homogenous medium with permittivities eQ,fio. According to this formula, 
we have P = eof^^^a^Eo, and M = ^J^^jga^Bo. With the condition e/eo,fio/fi < 1 for 
a hole in a perfect metal, we obtain the following two relations: P = — yO'^Eo, M = — a^Bg 



which are related by the same factor of 2 as for the Bethe case. 

The diffraction of light by a small aperture in a perfectly conducting plane can be regarded 
as the complementary case to the scattering of a wave by a small conducting particle. This 
point of view is in agreement with the electromagnetic Babinetjl7| theorem which yields this 
result directly ( see Smythe [oj. Landau 18|). Rayleigh has developed the theory of diffrac- 
tion by small particles and has found that the shape of the particle is not a fundamental 
parameter if one is interested in the far-field pattern of the scattering wave only. This is due 
to the fact that the light cannot distinguish the shape of the particle if the wavelength is very 
large compared to a. Therefore, we anticipate that the existence of two effective dipoles and 
their far-field behavior is very general and does not depend on the shape of the small hole or 
the screen geometries. The exact expressions of these dipoles may depend on the hole shape 
and their directions may differ from the Bethe case, but the factor of two between M and P 
as well as the presence of the fields and their dependance on the third power of the radius 
constitute probably a general result. A validity criterion for this statement can be found in 
the curvature of the surface at the hole location. If the curvature is too strong (i. e. if the 
curvature radius is too small), the hypothesis cannot be considered as reasonable due to the 
presence of strong boundary effects at the surface. Therefore, we expect the validity of the 
dipole model to diminish with the curvature radius. More precisely, the spatial domain in 
which the dipole model is applicable reduces progressively to a small solid angle around the 
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optical axis when the curvature radius decreases, i. e. when boundary effects invade space. 
We will see below that these general trends are confirmed by theory. 



III. THE QUASI-MULTIPOLE METHOD 

The method of quasi-multipoles developed here is the result of a convergence of three 
domains in electromagnetism. It lies at the intersection of i) the metallic wave guide theory 



for guides with variable cross sections 
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developed by Bouwkamp 
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2l|, ii) the spherical multipole expansion method 



|22|| for electromagnetic fields applicable to localized source 
distributions, and iii) the Hall 23| equilibrium solutions for conductors with conical geometry. 
For instance, the electromagnetic field radiated by a conical antenna or a conical hole can 
be described using this formalism. In order to simplify the discussion, we will use the same 
notations for spherical multipoles as those presented by Jackson |l6l|. Let us consider a 
source-free region of space with permittivities e and fi. A conical surface with half-angle P 
separates this space into an "oustide" region and an "inside" region. If 13 exceeds 7i/2, the 
external and the internal regions are mutually exchanged. For simplicity, these zones are 
considered as decoupled and totally independent. This is possible in the case of a perfectly 
conducting metallic surface in which tunneling of the field through the surface is forbidden. 
From now on, we will only consider the external problem. Maxwell's equations can be 
written (assuming an e"*"^* time dependence): 

E=^VxH ; B. = --^VxE 
ke kfx 

V.E = ; V.H = 0. (1) 

Making use of the identity (r-^) = r ■ (V^A) + 2V ■ A, these equations reduce to 

(V^ + k'^efi)r.E = ; (v^ + k\fi) r.H = 0. (2) 

By writing the Laplace operator in spherical coordinates, and separating the angular and 
radial variables, a typical solution of Eq. [2] can be expressed as follows: 

^ (r) = c^tmh^ {ky/ejir) Y^^^n {0, 4>) ■ (3) 

Both Hankel functions appear in this expansion ( the ± signs refer to outgoing and incom- 
ing waves, respectively). They satisfy the spherical Bessel equations for the radial variable 



r and the quasi-harmonic functions Yi^^m oc (cos 6*) ■ e*™''^ defined in the hterature [23l-l25|. 
The latter represent the generahzation to a conical geometry [2^ of the well-known spherical 
harmonics Yi^^ used in spherical eigenvalue problems. However, an important modification 
with respect to spherical harmonics is the presence of a non-integer parameter u which 
replaces the integer / and which depends on the boundary conditions on the cone. The 
possible values of u in Eq. [3] follow an infinite sequence which is a function of (3: Up = Vp 
{p integer). The exact relation depends on the boundary conditions on the cone. 
In analogy with the spherical case and with wave guide theory, we now decompose the 
electromagnetic field obeying a specific boundary condition on the cone into two parts: a 
magnetic quasi- multipole part (M) of the transverse electric field (TE), and an electric quasi- 
multipole part (E) of the transverse magnetic field (TM). These two parts are characterized 
by components of order (z/, m) satisfying 



-Y ^^.h^^\^r^)y.,rn{eA) (4) 



and 

r.i7j;^ = , r.Ej^„^ = (5) 

In this article we consider only outgoing radiation propagating in vacuum . For reasons 
of generality, we must use two types of indices ve {(3) and um (/?) because of a difference 
between the boundary conditions for TE and TM waves. For example, if the conductivity of 
the metallic cone is infinite (the perfect conductor), we must have the boundary conditions: 

n,,™ (^,0)|,__^ = O 



= 0, (6) 
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These expressions fix the authorized u values as a function of /3. Hence, an electromagnetic 
field satisfying specific boundary conditions can be written as an expansion in TE and TM 
waves, and we have 

H = V (kr) 

VE^rn ve \ ' / , "TT 



and 



E = V a*^ (kr) L^'-M.m (^^ 



E = ^ oa-Ea ; H = ^ aA-H;, (9) 

A A 

where A refers to TE and TM waves. In addition it can be demonstrated [l^ that we have 

ai.,-m = (-l)""^^ ■ (10) 

The coefficients a\ are determined by an integration with respect to the sohd angle which 
includes the entire allowed domain VLq 

ax = k'^ I r ■ (E X H^) rdVL. (11) 

Hence, we can express the angular distribution of radiated power, defined as ^ = ^'"^E x 
B*, as follows: 

dP_ ^ C / ■\i'i5+l JE) ^Y^E,m jd, (j)) 

■ " ^ ^ Jve {ve + 1) 



+ Hr'^' E ■ ipM^M= f. (12) 

i^M,m J Um [I'm + 1) 



In addition we can write the total radiated power as: 



^' Hhtf- (13) 



A 



IV. THE RADIATION PATTERN OF A COATED TAPERED FIBER TIP 



The tapered and metal-coated optical fiber tip used in NSOM|26l-l28| is characterized by 
a "funnel" geometry in which an aluminum layer (thickness of ~ 100 nm) is evaporated onto 
a glass cone truncated at its apex by a small disk (diameter of 2a ~ 60 nm, see Fig. |2]). 

The shape of the tip reveals a partially conical geometry. Therefore, the system is -to ffist 
approximation- equivalent to a conical antenna and it can be described by a quasi-multipole 



expansion. Nevertheless, in order to take into account the finite size of the aperture domain 
in our modehzation, we study the field outside a sphere of radius a/ sin /3 only. This so-called 
aperture zone represents the "terra incognita" near- field region (r/A < 1) of the fiber tip 
(see Fig. H]). Let e ~ yU 1 be the electric and magnetic permittivities of vacuum. The 
electromagnetic field in the region r /A > 1 can be described by an expansion using the quasi- 
multipole formalism. The coefficients a^^rn contain all properties of the electromagnetic field. 
In order to calculate E and B in the radiation zone, we need (see Eq. [TTj) to know the field 
on the spherical section VLq, i. e. in the near- field zone of the aperture. Unfortunately, there 
is as yet no complete theory of light diffraction by a small aperture in a screen with conical 
geometry, and we must rely on approximations. The ai^ ^n coefficients can be evaluated by 
using a Taylor series written in the vicinity of the tip origin. We obtain, using Eq. [TTl the 
following expressions to first order in r: 

ax ^ 27riP{-E (0) ■ — / r x (r x E^) r'^dn 

+ \ E[2A,a,/35/3^a (0) " Mx,a,pdpB^ (0)]. (14) 

In Eq. [HI the symmetric terms with prefactors Qa.o,/? and A^A,a,/3 are neglected in a first 
approximation. The integrals can be evaluated using the following properties: i) the /3 value 
is small and, as a consequence, the fraction (47r — VLq) / (47r) is negligible, ii) the authorized 
z/ values are close to 1, and iii) there are only the two first roots i'E,p=ii ^M,p=i associated 
with our boundary problem. These second and third hypothesis will be justified later on 
but they can be intuitively understood by realizing that in the limit of a very small tip angle 
/3, the theory must reduce to the case of a linear antenna in which the first and dominating 
radiation mode is the dipole term with / = 1. Hence, using Maxwell's equations, we have 

1 

af^,„ ^ -2«P(Eg,J.-^— ^E(O) 
atl^ ^ -2me{BtllJ.lj^,B (0) , (15) 
where the average is taken in the spherical domain of radius a/ sin/3. 

It can be noted that in the spherical harmonic expansion, the terms ||m|| > 1 are prohibited 
by recursion relations, but this is not the case here. 
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In order to calculate oa, we need to know the local field B (0) , E (0) at the center of the 
aperture. We would like to mention that a rigorous representation of the near-field can be 
drawn for a two-dimensional tip as shown in Fig. |31 Such a tip is described by magnetostatic 
and electrostatic potentials which can be found using the well-known complex potential 



method valid for the Laplace equation in two dimensions [29|. Numerical calculations by 
Novotny et al . [Soj] shows the same field behavior. The field topology in the 3D case can be 
anticipated from this 2D result. The choice of this field topology is based on the assumption 
that the mode entering the conical part of the tip is the one in the preceeding hollow, circular 
cylinder. This is principally the polarized TEll mode used in near-field microscopy. 
Turning back to the 3D case, we remark that, despite the lack of a rigorous analytical model, 
the near-field can be computed numerically using different methods. These calculations 
confirm the intuitive field topology obtained in our 2D result. Using our 2D model and 
numerical computations obtained by other authors, we postulate that the unknown near- 
field components are linked to each other by the "plane wave" relation : 

(E) = (B) X z, (16) 

where z is the cone symmetry axis. This condition may be intuitively understood by as- 
suming that the incident light propagating in the fiber is linearly polarized. This means 
that the symmetry of the field is conserved at the aperture as shown in Fig. [31 The elec- 
tric and magnetic amplitudes at the center of the aperture domain (i. e. at the top of 
the cone) can be determined from the Poynting theorem 18j which imposes the equality 

||B|p = ||E|p — / (E X B*) rrrffi in the spherical domain of the aperture zone. If 
the surface integral is equal to zero in this formula, i. e. if most of the energy is either 
transmitted or reflected but not stocked or dissipated in the near-field zone, we obtain 

(0) ^ S (0) as a good approximation. 
Hence, with this symmetry condition, the field is entirely determined and the calculation 
shows that the only non- vanishing terms have m = 1. The normalized angular distribution 
of radiated power is then 

2 

(17) 



1 dP ^2 
Pdn ~ 5 



Imi5simxf + 2Re>^-' 



'l^E (l^E +1) y l^M (l^M + 1) 

in which only the u values i>e and um have been considered . 

It is worth noting here that our derivation of Eq. [T7]can be obtained in a similar way by means 



of a quasi spherical eigenfunction expansion of the scalar Green's functions for a cone 3ll-l35| 



with Dirichlet (or Neumann type) boundary conditions depending on ve (respectively vm) 
and m. The scalar field -E (respectively r ■ B ), and consequently the coefficients ^ 
(respectively a*^^ ^ ) given by Eq. [TS]can be easily deduced from a surface integral on Qq. 
Conical Green's functions (which have already been used in the context of near-field optics 
for aperturless microscopes 

0,0) 

connect surface integrals, given by Eq. [11] and Eq. [151 



to the Huygens-Fresnel principle directly. 

In order to complete this section we can note than the ratio represents the physical 
quantity in measurements of the angular power emitted by actual fiber tips as explained 
later on. The total radiated power of the fiber tip can be expressed as a function of the field 
E(0) 



5c / ka 
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notal-5^^^j E(0r (18) 

This formula can be compared with the Bethe 18| result P = -^^^Eq, where Eq is the 



field at the center of the aperture. The same behavior appears in the two formulas: a 
dependence in and in E^. The difference arises from the (sin/3)^ which represents the 
second geometrical parameter of the fiber tip. It is important to note that this formula 
contains an unknown variable: the field Eq. This field is a function of the transmission of 
the tip which depends on the over-exponential decay in the conical wave guide due to the 



cut-off of the propagating mode in the fiber 38(] . This decay is not taken into account in our 



model which prevents us from computing the tip transmission. 

V. BOUNDARY CONDITIONS AND PERTURBATION METHOD 

In order to complete the solution of Maxwell's equations, we will describe the boundary 
conditions on the cone and justify the above assumption z/ ~ 1. A simple analogy with 
wave guide theory permits to establish the boundary conditions in a conical geometry (this 
corresponds to the case of a "wave guide" of variable cross section). For a perfect metal , 
we have 

Y.,,mie,<^)\, = 0, %^|, = (19) 

These equations impose a condition on ue and um which restricts the allowed values to 
infinite, growing sequences i'e,p and h'M,p that depend on (3 and m (p integer). Because 
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R = < 1 in the near-field, two successive roots u' < v associated with a TM mode (or 
a TE mode, respectively) obey: 



E,M 

E,M 
Qii/,m 



{kRY-''<^l. (20) 



This allows to neglect, to a good approximation, all roots except for the first two with p = 0. 
The calculation of these roots is possible for a small cone angle /3 using the approximation 
formula of Hobson and Schelkunoff 23|, |39| which give for m = 1: 

1 



1 . , 

Z^M,0,m=l ^ 1 + 2 , 6 • (21) 

For a tip angle of 15° we have uq ~ 1.03. This approximation can be compared with the 
values obtained by solving Eqs. [T9] numerically. Here we have for the same /3 ue^o — 1.033 
and UMfi — 0.967. 

In reality, aluminum is a good but not perfect metal which possesses a frequency-dependent 
complex dielectric constant ec = + ie". In the optical domain for A < 800 nm, the 
permittivity of aluminum can be described by the Lorentz formula: 

eA^)^l-^^- (22) 



The plasma frequency and the damping constant are given 3_8| by Up 15.56eV and 7 ~ 
0.608eV, respectively. At A ^ 633 nm, we have ec — — 54.4 + il8.8 which implies a refractive 
index ric = y/e = 1.18 + i7.12 and a skin depth 6 = ^ 7nm which is small but not 
negligible compared to the metallic clading of ~ 100 nm. The influence of the finite skin 
depth on the "perfect" and unperturbed TM and TE modes may be taken into account. 



The perturbed solution is discussed in the textbook of Jackson 



linking E and B on the surface with the aluminum permittivity 



16ILTI 
ty Q: 



The result is a condition 



E|, = E - ^ (e J) ~ ^ (23) 

This implies an effective current confined to the surface obtained by integration of Ohm's 
law in the aluminum layer of thickness 6 located just below the conductor surface. Using 
Eq. [231 we can establish 



1 kr dYy^^rn 
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BY 



de 



(24) 



As expected, this reduces to the unperturbed condition for a = oo. By using the substitution 

„ , (0) 

u — )■ z/*^") on the right-hand side of Eq. [2l] we have to lowest order ^Pe\s — Ie \^ ^^"^ 
^^l^ ~ /m'^m Is where the factors H/e.a/II are small. In order to calculate the perturbed 



dn 



V 



we use the two-dimensional Green's theorem in its integral form|l6[| and we obtain: 



i^f + l) - [ve + 1) ^ (1.4 + 20.23) kr 
(ym + l) - {vm + 1) =^ (1.4 + 20.23) 10"^^. (25) 

VEfi and vm,q are functions of two radii te — 18 nm and = A = 633 nm that are adjusted 
to reproduce the experimental results discussed below. The two values of z/ are complex and 
the real parts are close to unity as expected: ue — 0.95 — O.Oli, um — 0.96. These fits reveal 
that the adapted boundary conditions depend on a typical length r < A that is probably 
linked to the aperture radius a ^ 20nm (this value is in agreement with the aperture radius 
obtaiued b. sca.„,„g elect™ .nie.ose^yd)- W.th these values of . we can calculate the 



radiated power distribution of the tip|40l|. 



VI. COMPARISON WITH EXPERIMENTAL DATA 



Measurements of the far-field radiated power of tapered fiber tips in S and P plane 



have been carried out extensively by Obermiiller and Karrai ^, [s] . As seen in Fig. [6l 
there is an important backscattering effect in the P plane. The experimental data and the 
results obtained with the quasi-dipole model (using the values of i'e^m obtained above) are 



also shown in Fig. [6] for comparison (40l|. The agreement is good for most angles with the 
exception of extreme 6 values corresponding to grazing observation angles on the metallic 
cone. The disagreement between theory and experiment is probably due to higher terms in 
the expansion that we have neglected. This is confirmed by the recent independent results 
obtained by Shin et al.[41] which have fitted the experiment to an expansion in classical 
spherical multipoles. Their fit contains dipolar terms (where the factor 2 plays the same 
role as in our model) as well as quadrupolar (/ = 2) and octopolar {I = 3) terms. The 



behavior of the field in t 
due to surface roughness 



le extreme regions probably also depends on a modified diffusion 
42]. 
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The dipole model of Obermiiller and Karrai, which corresponds to z/ = 1 (a value close to 
the real values UEfi and I'Mfi), gives good results for the S polarization but cannot reproduce 
the backward emission in the P plane. The fact that the two values of i'e,m approach 
unity explains the efficiency of the dipole model for small azimuthal angles ( in this regime 
Y^ rri — ^i,m)- In addition we can formally justify the dipole model considering Eq. [TBI as 
the generalization for conical geometry of the fundamental definition of multipole moments 
which enter into the spherical harmonic multipole expansion, and which are used for localized 
source distributions. Using this analogy, we deduce the two effective dipoles associated with 
the aperture 

Meff = 2ixP,ff = 2ji=^B(0), (26) 

Once again, the characteristic factor of 2 between the magnetic and electric moments appears 
in this formula. Figs. H] and [5] are two equivalent representations of the spatial dependence 
of radiated power. In particular, the "orbital" representation of Fig. H] shows the important 
backward emission in the plane containing the effective electric dipole P and the tip axis 
(P plane). This backward effect is less important in the S plane containing the tip axis and 
the magnetic dipole due to the strong near-field effect on the metal surface. The factor of 
two in Eq. |26] (as in Eq. [T7|) stresses the difference between electric and magnetic dipoles 
as discussed in the introduction (namely an oscillating charge in the case of an electric 
dipole and a rotating charge for a magnetic dipole). This factor of two and the dependence 
of the dipoles on a and on Eq are similar to the Bethe model and the Rayleigh diffusion 
theory by subwavelength particles and holes, but the dipole orientations are different, a 
result unexplained by the Bethe-Bouwkamp solution. 

We can note in addition that the far fields radiated by a small distribution of electric and 
magnetic dipoles located in the apex zone which satisfy a relation of the type given by 
Eq. [261 are equivalent to those deduced from Eq. [15] (and consequently Eq. [TT]) in the limit 
of subwavelength distributions. This fact, which can be obtained by means of scalar Green's 
functions discussed above, gives further support to the analogy between dipoles and aperture 
at the apex of a cone. 

It is worthwhile to mention that A. Roberts 43| has used the method of quasi-multipoles 



for the interior problem of "small hole coupling of radiation into a near field probe". In 
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this problem the author used the hmit of small /3 where Y^^m (^j 0) ~ Ju[{y + |) 0]e^"^'^ and 
showed that the coupling of a wave into a conical wave guide implies effective dipoles given 
by the Bethe-Bouwkamp solution. These dipoles are perpendicular to the aperture, a result 
which is completely different from the problem of radiation by a tip which generates dipoles 
located in the aperture plane. 

For completeness, we wish to comment on some additional, interesting properties of the 
emission by conical tips that are relevant to the experiment. The quasi-multipole model 
shows that the aperture radius disappears in the formula Eq. [T7] and, consequently, in the 
radiation profile. We can understand this fact using the analogy with Rayleigh's and Mie's 



theory 12| of diffusion by small particles which follow the same behavior as our model and 
which can be associated with the impossibility to see details smaller than the wavelength 
of light. In particular, the aperture size cannot be a relevant variable in the limit of small 
apertures. For the same reason, the tip angle /3 does not play an important role in the 
diffraction profile. The dependence on the aperture radius is very different in the short- 
wavelength limit and for large apertures a > X. For a large hole, the transmission must 
approach unity and, by consequence, the aperture radius may appear in the profile . We 
have then ~ (fca)^ cos^ (6) (^^^^~^^)^) aiid P ~ cEqO^/S that can be compared with 
P ~ ^^1^2 Eq, valid for Bethe's theory and our conical subwavelength aperture model . 

VII. SUMMARY 

An optical fiber tip with a small aperture in the metal coating as used in NSOM is 
usually characterized by its angular transmission profile. In order to describe this profile 
theoretically, we have presented in this paper an analytical model for light diffraction through 
a subwavelength aperture located at the apex of a metallic screen with conical geometry. 
For this purpose, a field expansion in quasi spherical multipoles has been developed to solve 
Maxwell's equations analytically. Special care has been taken to consider realistic boundary 
conditions adapted not only to the conical geometry but also to the finite conductivity 
of the coating. Our model is able to reproduce and to explain the experimental results. 
The diffraction in the far-field is demonstrated to remain similar to that drawn with the 
2-dipole modelj^, [g] over a large portion of space. The latter model is formally justified in 
our approach. In particular, we find a characteristic factor of 2 between the expressions for 
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the magnetic and the electric dipole (which is reminescent of the Rayleigh diffusion theory 
by small particles and the Bethe problem). In addition, the previously unexplained large 
"backward" emission in the P plane (which does not exist in the S plane) results directly 
from our analysis. This "P emission" arises from a Poynting vector flow along the coated 
surface of the tip, an effect that cannot be described within a simple dipolar model. 
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All captions are in the order of appearence in the text. 




FIG. 1: A rigourous representation of the magnetic (a) and electric (b) fields in the vicinity of a 
circular aperture in a metallic screen. The aperture has subwavelength size and the screen has an 
infinite conductivity. The fields reduce to uniform far-fields Eq, Bq incident part of space (-), and 
to a dipolar field in the other part (+). 
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Far-field 




FIG. 2: An idealized metal-eoated fiber tip with its cliaracteristic funnel shape. The top angle (3 
is inferior to 15°, and the aperture diameter 2a is of the same order of magnitude as the aluminum 
coating thickness, i. e., 100 nm. 




(a) (b) 

FIG. 3: The electric (a) and magnetic (b) fields in the near-field zone of a two-dimensional cone, 
calculated rigourously with adapted complex potentials. The assumption is made that the topology 
of the 3D problem can be obtained by deforming the 2D case. In the 3D case, (a) and (b) planes 
are perpendicular to each other. In addition, the magnetic field in (a) and the electric field in (b) 
are perpendicular to both planes for symmetry reasons. 



18 



1,0 



(0 



0) 

c 

0) 



0,5 



0,0 



P-Polarization 



Quasi multipole 
model 




I I 



-180 -120 -60 60 120 180 

Angle (deg) 



1,0 



to 
c 

3 



3, 0,5 



CO 

c 

0) 



0,0 



S-Polarization 




Quasi multipoie 
model 



-180 -120 -60 60 

Angle (deg) 



120 180 



FIG. 4: "Orbital" representation of the angular radiated power of the tip {r = I {6,<p) ,6,<p ). The 
figure shows the important backward emission in the P plane passing through the tip axis and 
containing the effective electric dipole P. 
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FIG. 5: 3D isoline representation of the electromagnetic energy density (arb. units) in the far- 
field zone A;r S> 1. This figure shows the difference in the far-field energy distribution that exists 
between the S and the P polarization planes, and the presence of an important backscattering effect 
in the P plane. Additionally, the tip of the metallic cone is drawn in blue. The sphere at the tip 
apex corresponds to the (unknown) near-fied zone where the far-field divergent energy distribution 
oc ^ is no longer valid. 
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FIG. 6: Normalized angular distribution of radiated power for a) the P polarization where the 
incident electric field is parallel to the analysis plane, and b) for the S polarization where the 
detector is scanned perpendicularly to the plane of polarization of the incident light. Experimental 
data are shown in circles. The thick line corresponds to the present quasi- multipole model, the thin 
line to the 2-dipole model. The experimental data (aperture radius a = 20 nm, light wavelength 
A = 633 nm) are taken from Refs. 1,2. 
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